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Astigmatism of Skew Pencils in Optical Systems 
Containing Toric Surfaces 1 

Walter Weinstein 2 

A method is described for tracing skew rays through toric refracting surfaces and for 
determining the astigmatism of pencils around these rays. 

1. Introduction 

Toric refracting surfaces are used chiefly in two classes of optical systems, spectacle lenses, 
and anamorphotic photographic and projection objectives. In spectacles the toric surfaces 
are arranged to correct the axial astigmatism of the eye, and in anamorphotic systems two or 
more toric (possibly cylindrical) surfaces are included to give stigmatic imagery with different 
magnifications in different directions in the image plane. 

When the axial pencil is refracted by a toric, 3 some astigmatism is introduced; if a pencil 
is incident obliquely, either centrally or eccentrically, astigmatism of different magnitude is 
introduced. It is therefore desirable in the case of both types of optical systems mentioned 
above to determine the astigmatism of obliquely incident pencils, in order to assess the imagery 
in the outer parts of the field. This is the problem which was studied in the work to be described. 

The problem arose in connection with spectacle lenses and may be outlined as follows: 

Figure 1 shows front and side views of a spectacle lens with toric surface ; is the center of 
the lens, C is the center of rotation of the eye, and P u P 2 , P 3 are points of incidence of the chief 
rays of oblique pencils. Let the meridian and equatorial sections of the toric be, respectively, 
in the vertical and horizontal planes 0P X and 0P 3 . Then the problem for chief rays passing 
through the axis and through P x and P 3 is trivial. The chief ray itself can be traced by the 
ordinary ray-tracing formulas for meridian rays and spherical surfaces, and the astigmatism 
can be determined by the usual methods for close sagittal and tangential fans, with the slight 
difference that the curvatures of the refracting surface in the two directions are not equal 
(see eq (8)). 




Figure 1. 



If, however, the point of incidence of the chief ray does not lie in the meridian or equatorial 
section, as, for example, P 2 in figure 1, the problem becomes more complicated. In the first 
place, even if the incident chief ray intersects the axis of the system, the refracted ray in general 
does not, so that the chief ray must for complete generality be considered to be skew; thus the 
first part of the problem is to trace a skew chief ray through the toric. 

1 The work described in this paper was carried out while the author was a guest worker in the Optical Instruments Section, National Bureau 
of Standards, from January to July 1954. 

"■ Present address: Imperial College, London, England. 

3 In what follows, "toric surface" will be abbreviated to "toric." 
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Second, with oblique incidence as at any of the points P u P 2 , or P 3 , astigmatism is produced 
in the refracted pencil for two reasons: (a) because the angle of incidence is not zero, as in the 
ordinary case of spherical surfaces, and (b) because the refracting surface has two different 
principal curvatures at the point of incidence. In the case of P 2 the principal directions of the 
astigmatism due to (a) are parallel and normal to the plane of incidence of the chief ray, i. e., 
approximately parallel and perpendicular to line, 0P 2 ; the principal directions of the astigma- 
tism due to (b) are, however, approximately parallel to 0P\ and 0P 3 . Thus the resultant 
astigmatic pencil will have its focal lines in some intermediate directions, and the determina- 
tion of these directions further complicates the calculation. 4 It follows that the formal dis- 
tinction between sagittal and tangential directions can no longer be made; however, for the 
configurations of chief ray and toric that occur in spectacle lenses the rotation of the focal lines 
is unlikely to be above 10°, so that there is little difficulty in relating the results for skew chief 
rays to those for rays in either principal section. 

The results of such calculations would constitute a survey of the astigmatism of a toric 
spectacle lens at all obliquities and azimuths. They would help to answer questions such as the 
following: What is the best bending of a toric spectacle lens? What range of powers can be 
manufactured with a common back curvature without departing appreciably from the best 
bending? Is there any advantage in making both surfaces of a spectacle lens toric in shape? 
Such questions as these arise in the design of so-called best-form spectacle lenses. These have 
previously been designed on the basis of computations in the two principal sections only. 5 6 
However, an illustration is given by Henker 6 showing the quality of an image produced by a 
pencil outside the principal sections of a typical sphero-toric lens of the Punktal best-form 
series; this illustration suggests that there is still room for improvement in the image. The 
work to be described makes it possible to compute the imagery by skew pencils and so determine 
whether any improvements can in fact be made. 

Section 2 is a recapitulation of the properties of astigmatic pencils that are required in 
deriving the astigmatism equations. In section 3 the equations of a toric are given, and other 
properties are derived, in section 4 the ray-tracing equations for the skew chief ray are derived, 
and in section 5 the astigmatism formulas are derived. Section 6 deals with the forms in which 
the input data must be given. In section 7 the mode of presentation of the results is discussed, 
the results of a numerical example are given, and the possible uses of the computation are 
discussed. 

2. Astigmatic Pencils 

For the purposes of the present paper a pencil will be considered to be made up of a narrow 
bundle of rays originating from a point in the object and surrounding a chief, or principal, ray, 
together with all the corresponding wave fronts. The wave fronts are surfaces of constant 
optical path length from the object point, and the rays are normals to the wave fronts. 

With a rectangular coordinate system with the x axis along the principal ray, the equation 
of the wave front passing through the origin will take the form 

x=c 2 oy 2J rC n yz J rC m z 2 + terms of higher degree in y and z; (1) 

linear terms are absent because the yz plane must be tangent to the wave front at the origin. 
Also it is possible to rotate the y and z axes about the x axis, so that the coefficient of the term 
in yz vanishes when the equation takes the form 

x=i {c r y 2j rc s z 2 ) + higher-degree terms. (2) 



4 This difficulty does not arise for points on the principal sections of the toric, such as P- and Pz, because then the directions of the focal lines 
are known by symmetry to be parallel and perpendicular to the principal sections. 

5 H. H. Emsley and W. Swaine, Ophthalmic lenses, 6th ed. (London, 1951). 

6 0. Henker, Introduction to the theory of spectacles (Translated by R. Kanthack) (Jena, 1924). 
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Effects depending on terms in y and z of higher degree than the second will be neglected. 
To this degree of approximation it is easily shown that all rays of the pencil; i. e., normals to 
the surface (2) intersect one line through the point (l/c TJ 0, 0) parallel to the z axis and another 
line through the point (l/c„ 0, 0) parallel to the y axis, these being the two astigmatic focal 
lines. Thus this approximation takes account of astigmatism but not of aberrations such as 
spherical aberration and coma, which depend on higher powers of the aperture of the pencil. 
The coefficients have the significance that c r is the curvature of the section of the wave front by 
the xy plane and c s , that by the xz plane 

Some results from differential geometry are now required. At any nonsingular point P 
of a surface the normal curvature (i. e., the curvature of the section of the surface by a plane 
through the normal) has a maximum and a minimum for two section planes which are at right 
angles to each other; furthermore, if these two planes, together with the tangent plane at P, are 
taken as coordinate planes, the equation of the surface in the neighborhood of P takes the form 
of eq (2), that is, the yz term is absent. The maximum and minimum normal curvatures are 
called the principal curvatures, and the corresponding section planes meet the tangent plane 
at P in lines that lie along the principal directions at P. Thus the principal directions are along 
the y and z axes for eq (2). 

It is also known, by Eider's theorem, that if a normal section plane makes an angle a with 
the xy plane, the curvature in this section is 



c(a) =c r cos 2 a-\-c s sin 2 a. 



(3) 



From this it follows that these quantities are required to characterize an astigmatic 
pencil: 

(a) The direction of the chief ray; this will be given by a unit vector, Q, with components 
X, Y, and Z. 

(b) The principal curvatures, c r and c s , of the wave front at a point on the chief ray. 

(c) The principal directions on the wave front; these will be given by unit vectors R and S, 
and from the above it follows that Q, R, and S are all perpendicular to each other. 

3. Toric Surface 

A toric surface is generated by the rotation of a circle (the generating circle, abbreviated 
to G. C.) about an axis in its plane but not passing through its center, as in figure 2, a and 
b. In optical applications the lens surface is, of course, only a cap-shaped portion of the 
whole surface, roughly centered on the point A. Such a surface is alined with a centered 
system of spherical surfaces by making the optical axis lie in the normal plane to the rotation 
axis, which contains the center of the G. C. The optical axis then meets the toric in the 
pole A. 



ROTATION AXIS 



ROTATION AXIS 




RING TORIC 



Figure 2. 




BARREL TORIC 
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The toric refracting surface has different powers (in the optical sense) in different sections- 
through the optical axis: the powers are determined by the radius of the G. C. for the section 
in the plane of the rotation axis and by the distance from A to the rotation axis for that 
perpendicular to it; these are known as the meridian and equatorial sections, respectively. 

It can be seen that for given powers in the principal sections there are two possible forms 
of toric surface: the radius of the G. C. can be either less than or greater than the distance 
from A to the rotation axis, the two forms being known, respectively, as ring- and barrel- 
shaped torics. The ring form seems to be used almost exclusively for spectacles in commercial 
practice, but to the best of the author's knowledge, this is merely because more rings than 
barrels of given curvatures can be surfaced in one operation, not for any optical advantage. 
The possibility of improvements from the use of the barrel shape seems not to be considered 
in the literature on "best-form" series. 

3.1. Equations of the Toric 

One can write down the equation of the toric surface by taking rectangular coordinate 
axes with the origin at the pole A of the toric refracting surface (fig. 3), the x axis along the 
optical axis, the y axis parallel to the rotation axis, and the z axis in the equatorial plane. 
Let r y and r z be the radii of curvature in the xy and xz sections, respective^. Take angular 
parameters, 6 and <£, as in the figure, so that 6 is the longitude and </> the latitude of a point 
on the toric. Then the parametric equations of the toric are 



x=r z — (r z — r y -\-r y cos </>) cos 6 
y=r y sin 4> 

z=(r z — r y J rry cos </>) sin 6 



(4) 



It can be seen that eq (4) represent a ring- or a barrel-shaped toric according as r y is less 
than or greater than r z in absolute magnitude, 7 in both cases the rotation axis is parallel to the 
y axis. When the rotation axis is to be parallel to the z axis the equations are obtained be 
interchanging y and z throughout, including subscripts, in eq (4); this is true also for the 
equations in the remainder of section 3. 
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An explicit equation can be obtained by eliminating 6 and <£ in eq (4) and solving for x: 

x=r z ± { [r z -r v ±{r y *-yy*Y-z 2 }*. (5) 

The ambiguous signs are due to the fact that there are four (possibly complex) values of 
x for any given y and z, since the toric is a quartic surface; the required solution is that which 
tends to zero as y and z tend to zero, so that the correct signs are, respectively, minus and 
plus. When these are inserted the result is unsuitable for numerical computation when y and 

7 Provided r y and r t have the same sign, which is generally the case in spectacle lenses. 
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z are small compared to r v and r z , as is generally 1 Ik 4 case. This difficulty can be circumvented 
by transforming- the equation by means of conjugate surds in the usual way and using curva- 
tures instead of radii, giving 

F 

x= G{G+[G 2 -c z F]v>Y 
where 

F=c 2 z*G 2 + 2c y tfG-cic z y\ G= 1 + (1 -c v Y) * 

In this form the equation of the toric is suitable for numerical computation with all 
relevant ranges of the variables. 

3.2. Principal Directions and Curvatures 

In considering the refraction of an astigmatic pencil through a toric, it is necessary to know 
the shape of the toric in the neighborhood of the point of incidence, Thus, as in section 2, the 
following are required: 

(a) The direction of the normal at the point of incidence; this will be given by a unit 
vector, U, with components, X, /x, and v. 

(b) The principal curvatures, c v and c w , at the point of incidence. 

(c) The principal directions at the point of incidence; these will be specified by unit vectors, 
I 7 and W, and again the vectors, U, V, and W, form a trirectangular system. These quantities 
are easily determined for a current point (x,y,z); if the toric has principal curvatures, c y and c Z} 
at the origin and if the rotation axis is parallel to the y axis, as in section 2.1, the results are 

U= (cos <t> cos 6, — sin </>, — cos sin 0) , (7) 

Q • Gc z z 

<£ = arcsin c y y, 0=arcsm -^ — — 2 > 



-CzCylf 



C v Cy 



_Gcz cos 4> Y (8) 

Cw -G-c^f J 



V= (sin </> cos 6, cos </>, — sin sin i 

(9) 
W= (sin 6, 0, cos 6) J 
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4. Tracing the Skew Chief Ray 

Using the quantities defined in sections 2 and 3.2, the problem can be stated as follows: 
Given Q and the coordinates of a point on the chief ray, the coordinates (x,y,z) of the 
point of incidence of the ray on the toric must be determined (transfer process); next U at the 
point of incidence must be found and Q '\ the direction of the refracted chief ray, determined 
(refraction process). The quantities c r , c s , R, and S will be known for the incident pencil 
and c Vy c W) V, and W can be found, since the point of incidence is known; from these the 
quantities c' r , c' s , R', and S' for the refracted pencil must be determined. This last, the astig- 
matism calculation, will be described in section 5. 

4.1. Transfer Process 

The transfer problem is that of finding the intersection of a straight line with a toric. A 
direct approach would require the solution of an algebraic equation of the fourth degree; this 
would be a laborious process, so an iterative method was adopted. In addition to the general 
advantage of iterative methods, that mistakes in the early stages of the calculation are auto- 
ma tica II y eliminated at the cost of extra computing time, this method also selects the correct 
solution from the four possible solutions. 
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Let Pi be a point on the ray near the required point of incidence. 8 A point P[ on the 
toric is found with the same y and z coordinates and the tangent plane T x at P[ is determined. 
The next approximation P 2 is taken at the intersection of the ray with 7\ and the 
cycle P 2 -^P2-^T 2 ->Pz . . . repeated until the coordinates of P n and P n +i are the same to the 
required number of significant figures. These are then the coordinates of P, the point of 
incidence. 

Let the coordinates of P n be (x ny y n ,z n ) and those of P' n be (x' n ,y n ,z n ); then x n is obtained 
by substituting y n and z n in eq (6). The equation of the plane T n is K(x—x n ) + ii n (y—y n ) + 
v n (z—z n ) = 0, 
where \, /x n , and vz are given by eq (7) so that x n+1 is given by 



x n + 1 — x n -+- b n 

\ n X{X n X n 



d n = 



(10) 



Also, 



\ n X+^ n Y+v n Z 



Y "} 

y n+1 =y n +j^d n 



(11) 



This completes the iterative cycle. It can be shown to be a second-order iteration, 9 that is, 
5 n+1 is of order of magnitude 8 2 n , so that it is always rapidly convergent, provided the initial 
approximation is good enough. 

The choice of the initial approximation depends on the available means of computation. 
If a high-speed automatic-computing engine is available, Pi may be taken as the point where 
the ray meets the yz plane. Let (x ,y ,z ) be the coordinates of the point P where the ray 
leaves the previous refracting surface (or the entrance pupil); 10 then Pj is given by 



#i=0, 

yi=yo-x Y/X, 
Zi=zq—x Z/X j 



(12) 



Starting from this approximation, some four or five iterations would probably be needed to 
reach a result correct to six significant figures. 

If a desk computing machine is to be used, an appreciable saving in computing time can be 
effected by using an improved first estimate. It usually happens that chief rays directed at 
the same pupil have already been traced in the meridian and equatorial sections of the toric. 
Let the x coordinates of their points of incidence be x y and x z , respectively. The section of 
the toric by a plane through the x axis and making an angle a with the xy plane is approximately 
a circle, of curvature c^cos 2 a-{-c x sin 2 a, so that a chief ray in this plane will meet the toric at a 
point of which the x coordinate is approximately x y cos 2 a-\-x z sin 2 a. This Pi is given by 



Xi=x y cos 2 OL+Xz sin 2 a 

yi=yo+(xi—x )YIX 

zi=z J r(xi—xo)ZIX 



(13) 



This approximation is usually so good that only one iteration is needed to obtain six 
significant figures in the coordinates of P. 

8 The method of obtaining Pi will be discussed below. 

9 W. Weinstein, Proc. Phys. Soc. (London) [B] 65, 731 (1952). 

10 (jo, yo, zo) are taken in the coordinate system referring to the current surface, not that on which Po lies. 



256 



4.2. Refraction Process 

The refraction calculation for the chief ray can now be carried out by well-known formulas. 
Let / be the angle of incidence and n the index of refraction of the medium in which the ray is 
incident; let primes denote corresponding refracted quantities. Then 



cos/-Q.f/, 



cos r 



--+H 



(1 — cos : 



\l/2 

i)) . 

Q'=~ Q+fcos I'-—, cos i) U 



(14) 
(15) 

(16) 



In eq (14) QJJ denotes the scalar product of Q and V. 

This completes the refraction process and the coordinates and direction cosines of the 
refracted ray are then used, with a suitable change of origin, for the incident ray at the next 
surface. 

5. Astigmatism Calculation 

In order to determine equations giving R f , S', c' r , and c' s , in terms of known quantities, a 
general theorem is required. Referring to figure 4, let P be the point of incidence of the chief 
ray, and let B be the point of incidence of another ray of the pencil. Suppose the incident and 
refracted wavefronts through P to be constructed, and let the ray through B meet these in C 
and C'j respectively. 11 Then, since the optical path length between two wavefronts is the same 
along any ray, it follows that nGB=n'C'B' or, if CB is denoted by q, 

A(:nq) = 0. 

Here the symbol, A, means that the refracted quantity minus the incident quantity is to be 
taken. This relationship is true for any relative positions of P and B and for any shape of 
refracting surface. 

In order to use eq (17) for the present purpose, q and q' will be expressed in terms of the 
coordinates of the point B on the refracting surface and the parameters of the wavefronts. 
The resulting equation will give the normal curvature of the refracted wavefront in a section 
containing the ray through B and by differentiating with respect to the azimuth of B the maxi- 
mum and minimum normal curvatures, c' r and c' s , can be found. The same calculation gives 
the directions R' and S' of the principal sections. 



INCIDENT REFRACTED 

WAVEFRONT WAVEFRONT 




CHIEF RAY 



Figure 4. 



The part of the refracted wavefront shown in the diagram is virtual. 
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5.1. Normal Curvature of the Refracted Wavefront 

Figure 5 shows the refracting surface near the point of incidence. Let P and B have the 
same significance as in figure 4 and let Fand Wbe the unit vectors along the principal directions. 
Let Bi be the projection of B on the VW plane; let a be the angle between PB X and V, and 
let PB be denoted by p. Then the position of B relative to P will be given by p and a as 
coordinates. 




Figure 5. 

In order to express q in terms of p and a, use is made of the fact that CB in figure 4 is a 
normal to the wavefront. It is found (see section 5.1, a) that 

q=p(Q-V cos a + Q-W sin a) + ^ p 2 {Q-U(c v cos 2 a+c w sin 2 a)— c r (R-V cos a + R- W sin a) 2 

-c s (S-V cos a+S-W sin a) 2 } + terms in p\ (18) 

and there is a similar equation for q\ 
Substituting for q and q' in eq (17), 

pAn(Q-V cos a+Q-W sin a) + lp 2 An{Q-U(c v cos 2 a+c w sin 2 a)— c r (R'V cos a+R-W sin a) 2 

— c s (S-V cos a+S-W sin a) 2 } + terms in p z =0. (19) 

Only terms up to and including the square of p are taken since, as explained in section 2, 
this degree of approximation takes account of astigmatism. 

Equation 19 holds good for any values of p and a, so that the coefficients of different 
powers of p may be separately equated to zero. If this is done, the linear term yields an altern- 
ative form of SnelFs law and on putting a=0 and tt/2, the equations 

A[nQ-V]=A[nQ-W] = Q (20) 

are obtained. These could serve as a check on the refraction calculation. 

Equating the coefficient of p 2 in eq (19) to zero and dividing through by 1 — (Q' V cos a + 
Q'-W sin a) 2 gives 



cffi 'V cos a+R'.W sin a) 2 +c' s J§'-V cos a+S'-Wsina) 2 

(21) 



l-(Q'-Fcos a+Q'-Wsm. a) 2 
nc T {R-V cos a+R- W sin a) 2 +nc s (S- V cos a+S- W sin a) 2 +(c v cos 2 a-\-c w sin a)AnQ-U 



n'[l-(Q'-Vcos a+Q'-W sin a) 2 ] 

Now the equation of the refracted wavefront can be shown to be 

[e'r(R'-V cos a+R'-W sin a) 2 +c' s (S'-V cos a+S'-W sin a) 2 } (x?+tf) 
x (R'-Vcos a+R' W sin a) 2 +(5'Tcos a+S'-JFsin a) 2 

taking the x, y, and z axes along Q', ft', and S'. Thus the left hand side of eq (21) is c'(a), the 
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normal curvature of the refracted wavefront in a section containing the point B with angular 
coordinate a on the toric. 

a. Derivation of Equation (18) 

Two auxiliary coordinate systems are required, each with origin at P, the point of inci- 
dence. One, the x, y, z system, has its axes respectively parallel to Q, R, and S; the other, 
x', y', and z f has its axes parallel to U, V, and W. The two are related by the equations 

x=(Q-U)z'+(Q.V)y' + (Q.W)z', y=(R-U)x'+(R.V)y'+(R.W)z', 

z=(S.U)x'+(S.V)y' + (S-W)z\ 

If the point B on the refracting surface (fig. 4) has polar coordinates (p, a), its coordinates 
in the x', y', z' system are 12 

x' = i(c v cos 2 a-\-c w sin 2 a)p 2 +0(p*), y'=p cos a-{-0(p s ), z'=p sin a + Oip 3 ). 

Let B have coordinates (u, v, w) in the (x, y, z) system; then 

u =%(Q.U)(c v cos 2 a + c w sin 2 a)p 2 +(Q-V)p cos a+{Q-W)p sin a + 0(p 3 ), 

V=(R-V)p cos a+(R-W)p sin a + 0(p 2 ), 

w= (S> V)p cos a+ (S- W)p sin a+0(p 2 ). 

Let (x, y, z) be the coordinates of C in the x, y, z system, so that 

x=h(c r y 2 +c s z 2 ) + 0(y\z*). 

The direction cosines of the normal to the wavefront at C are, to the approximation 0(y 2 ,z 2 ), 
(l, — c r y,—c s z), and since CB is a ray, or wavefront normal, it follows that 

u=x~\- q, v=y — c r yq, w=z—c s zq. 

There are now enough relations to eliminate x, y, z, u, v, and w. The resulting equation involv- 
ing p, q, and a can be solved for q, remembering that terms which are 0{p i ,(f) can be neglected, 
since only astigmatism is in question. The solution for q is found to be eq (18). 

5.2. The Principal Curvatures of the Refracted Wavefront 

The principal curvatures c' s and c' r are the maximum and minimum values of c'(a). 
Since all quantities on the right of eq (21) are known with the exception of a, c' r and c' s can be 
obtained in the usual way by differentiating with respect to a and equating to zero. Let 

a=±[nc r ((R.Vy+(R-Wy 2 ) + nc s ttS.Vy+ (22) 

b=\[nc T {(R-V) 2 -{R-W) 2 )+ncX(S-V) 2 -(S^W) 2 )+ (23) 

c=nc r (R- V) (R W)+nc s (S- V) (S- W), (24) 

d=K(?-(Q''Vy-(Q'-W) 2 ), (25) 

e= 1 2«Q'-V) 2 -(Q'.W) 2 ), (26) 

f=-(Q'-V)(Q'-W). (27) 

The equation for a is 

(ea — bd) sin 2 a+idc — af) cos 2 a=(bf—ec), 



'-' The expression 0(p n ) used in this section denotes terms of order of magnitude p n . 
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and the two solutions are 



where 



7-/3 7T 7+/ 



and 



iS=arcsin 



7=arcsin 



dc—af 



+ [(ea-bd) 2 +(dc-af) 2 ]^ 

bj—ec 

[{ea-bd) 2 +(dc-aj) 2 \^ 



(28) 
(29) 
(30) 



The values of c[ and c' r are then obtained by putting a = a x and a 2 in the right hand side of 
eq (21). 

5.3. Principal Directions on the Refracted Wavefront 

In figure 6, let P be the point of incidence, PB the refracting surface, and PC the refracted 
wavefront. Let Q' be the unit vector along the chief ray and let C'B be another ray of the 
astigmatic pencil. Let the projections of C and B on the tangent planes to the wavefront and 
refracting surface be (\ and B x respectively. Then V cos a+ W sin a is a vector in the direction 
PB X , so that a unit vector in the direction PC/ is given by 



F cos a+W sin a-Q'(Q'.V cos a+Q'-W sin a) 

~ [1-(Q'.V cos a+Q'-W zm a) 2 ] 1 ' 2 



(31) 



Thus, the vectors R' and S' are obtained by putting a=oti and a 2 in the expression (31). It 
should be noticed that in this derivation it was assumed that the points C and B lie on a plane 
through Q f ', i. e., that the ray through C intersects the chief ray. This is not strictly true in 
general, but the error depends on terms in the equation of the wavefront of the third degree in 
the aperture, so that astigmatism calculations are not affected. 



Vcos a tW slna 




Figure 6. 



This completes the refraction computation for astigmatism. The R' and S' vectors are 
transferred without change to be the new R and S vectors at the next surface. In order to 
transfer the curvatures the distance L between surfaces along the chief ray is required; this 
is easily found to be (x—x Q )/X, in the notation of section 4.1. Then the new curvatures are 
given by 



(l/c f ) +1 =(l/c r )'-2^ 



(l/C,) +I = (l/c s )'-ij 



y 



(32) 



5 Of the possible solutions for /9, one must be chosen which satisfies 



/8=arccos~ 



ea—bd 



+[f>a-6rf)'+(dc-a/)2] 



-+a multiple of 2w. 
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6. Input Data 

If the wavefront incident on the first surface is not already astigmatic, some arbitrary 
directions for R and S must be assigned. Similarly if a refracting surface is spherical, direc- 
tions for V and JFmust be assigned. 

These pairs of vectors must form trirect angular sets with Q and U, respectively, but subject 
to this condition they can be chosen to suit computational convenience. Choosing R to be 
coplanar with Q and the y axis, it is found that 



R 



S= 



(-xr,i-r 2 ,-rz) 

(1-F 2 ) 1/2 
(Z,0,-Z) 



(1-F 2 ) !/ 



(33) 



Since Y 2 is rarely greater than 0.1, the vectors so defined are well determined numericall}-. 
In a similar way V and W can be taken as 



V= 

w 



(— X^,l— m 2 , — M^ / ) > 



(1~M 2 ) 1/2 

OsQ,-\) 



(1-M 2 ) 1/2 



(34) 



It is also necessary to choose between ring and barrel shapes for the toric surfaces. In the 
equations of section 3.1 the rotation axis of the toric is taken parallel to the y axis and the toric 
will be ring or barrel shape according as c y is numerically geater or less than c z . Thus, if it is 
required to have a barrel shape with |c y |>|c 2 | or a ring shape with K|<C|c 2 |, the equations 
must be changed so that the rotation axis is parallel to the z axis. This can be done by inter- 
changing y and z in the symbols and subscripts of eq (4) to (8), inclusive, as mentioned in 
section 3.1. 

7. Presentation of Results 

7.1. A Numerical Example 

In the application to spectacle lenses the results are presented in the form of a diagram 
giving the astigmatism of a number of pencils and the directions of the principal sections on 
the wave front. An example is given in figure 7. 



5.730 

3509 




5.730 
3.624 




RING TORIC BARREL TORIC 

Figure 7. 
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The wavefront curvatures are marked on the appropriate sections, in diopters with respect 
to the vertex sphere. The apparent rotation of the principal sections of the wavefront can be 
found from the consideration that the vertical plane containing the chief ray contains the 
direction which is apparently vertical to the eye. This plane is the plane through Q and 
parallel to the y axis. From this it can be shown that the angle, A, through which the R and 
S vectors appear to be turned is given by 

tan^l=-, (35) 

s 2 

where r 2 and s 2 are the y components of R and S. 

The nominal powers of the lens in this example were 5.5 diopters in the vertical section 
and 3.62 diopters in the horizontal section. The curvatures of the front (toric) surface were 
+ 0.2199 cm -1 and +0.1839 cm" 1 , that of the back (spherical) surface was —0.1147 cm" 1 , the 
center thickness was 0.48 cm and the refractive index was taken as 1.523. The chief ray was 
made to intersect the axis 27 mm behind the lens at an angle of 30°. The results for both ring 
and barrel shapes with the same nominal curvatures are given. 

Computations have also been done on the prescriptions of a number of lenses and measure- 
ments have been made on these with a vertex power measuring instrument. The agreement 
is satisfactory within the limits of experimental error. 

7.2. Discussion 

In the numerical example given above the rotation of the focal lines in intermediate 
azimuths is negligible; this may be expected to be the case for all but very heavily cylindrical 
spectacle lenses of best-form type. It can also be seen that the variation of astigmatism over 
the field is less fo^ the barrel than for the ring shape (0.14 D and 0.21 D, respectively). This, 
of course, is not a general conclusion, but it suggests that further computation would be 
desirable. It is possible that for heavy cylinders a considerable improvement might be 
obtained by using barrel shapes, but this might be outweighed by the increased manufacturing 
costs. This would depend also on visual tolerances for astigmatism; it seems that more 
experimental work is required on this subject. 

A different question is that of the accuracy with which torics are produced. The principal 
curvatures check within the limits of accuracy of a lens gage, and the sagittal curvatures change 
or remain constant along the meridan or equatorial sections, respectively, to the same accuracy. 
However, this may not be enough to prove that the lens surface is a true toric to an accuracy 
which justifies the calculations for chief rays in intermediate azimuths. The mode of manu- 
facture is certainly such as to produce a very close approximation to a toric if the surfacing 
machine really acts in the intended manner, but it is quite possible that distortions and non- 
uniform wear in the tools and machine parts may cause departures from the toric shape. No 
measurements other than checks on the principal curvatures seem to have been made. A 
detailed check on the x coordinates of a network of points would be needed. 



Thanks are due to my friends and colleagues of the Optical Instruments Section for their 
interest in this work, in particular O. N. Stavroudis and Loyd Sutton; the latter did the 
numerical computations and checked the mathematics. 

Washington, June 29, 1954. 
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